Discontinuous Electromagnetic Fields Using Huygens Sources For Wavefront 

Manipulation 
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We introduce the idea of discontinuous electric and magnetic fields at a boundary to design and 
shape wavefronts in an arbitrary manner. To create this discontinuity in the field we use electric and 
magnetic currents which act like a Huygens source to radiate the desired wavefront. These currents 
can be synthesized either by an array of electric and magnetic dipoles or by a combined impedance 
and admittance surface. A dipole array is an active implementation to impose discontinuous fields 
while the impedance/admittance surface acts as a passive one. We then expand on our previous 
work showing how electric and magnetic dipole arrays can be used to cloak an object demonstrating 
two novel cloaking schemes. We also show how to arbitrarily refract a beam using a set of impedance 
and admittance surfaces. Refraction using the idea of discontinuous fields is shown to be a more 
general case of refraction using phase discontinuities. 

PACS numbers: 42.25.Bs, 41.20.-q 



INTRODUCTION 



One of the most basic boundary conditions in elec- 
tromagnetic theory concerns the interface between two 
different materials. As propagating waves traverse the 
boundary they are reflected and refracted with different 
amplitudes. This occurs specifically because the tangen- 
tial electric and magnetic fields across the boundary are 
conserved. One way to interpret this physical behaviour 
is that by maintaining the continuity of the tangential 
fields across a material boundary one can subsequently 
alter and control the propagation of electromagnetic 
fields in a region. This fact forms the basis for 
engineering and designing many electromagnetic and op- 
tical devices, including waveguides, lenses and scatterers. 

Because this idea is so fundamental to electromag- 
netic theory, the continuity of the fields across interfaces 
is often taken for granted. However it is worth asking 
a simple question. What if instead of using material 
interfaces, and thus continuous fields, to engineer 
desired wave propagation effects we could somehow 
impose discontinuities in the field, and thereby alter the 
propagation and amplitude of the wave. In this scenario, 
a desired discontinuity in the electromagnetic field is 
introduced across a boundary by purposely engineering 
the boundary itself. Here, we are building on work done 
in [1] where we introduced a discontinuous field at a 
boundary to cloak an object. We now seek to generalize 
this idea by exploring in more detail how such bound- 
aries are realized, and how can they be used to alter and 
control the electromagnetic field in a variety of scenarios. 

Specifically, we approach this problem by propos- 
ing a general way to impose a discontinuity using 
currents at an interface. These currents will act like a 
Huygens source which radiates the desired wavefront to 



create the discontinuity in the electromagnetic field. We 
first propose two ways to implement our discontinuity, 
one active and one passive. Using the active method 
we synthesize some novel cloaking examples building on 
our previously developed cloaking approach in [1]. With 
the passive method we develop a novel way to refract 
a plane- wave. Throughout these specific examples we 
will demonstrate how discontinuous fields can be used 
to create novel electromagnetic devices. 



DISCONTINUOUS ELECTROMAGNETIC 
FIELDS 




FIG. 1. An arbitrary boundary in free space upon which a 
discontinuity exists between the electric and magnetic fields 
on either side of the boundary. Electric and magnetic currents 
are imposed on the boundary to support the discontinuity. 
This is the motivating idea behind altering the fields using 
electric and magnetic currents. 

The answer to the question of how one imposes 
a discontinuity in the electromagnetic field, follows 
from a simple examination of the boundary conditions 
presented in electromagnetic theory. To see this let 
us imagine an empty space (free space) without any 
materials present. We then have an arbitrary boundary 
that divides our space into two halves and is described by 
a unit normal n as shown in Fig. 1. The two half spaces 



are also numbered respectively as shown. Without loss 
of generality we assume (throughout this paper) a field 
distribution that is invariant along the z-axis (J^ =0). 
We also assume a polarization of the fields with the 
electric field along the z-axis, E = Ez, and the magnetic 
field in the plane. 




FIG. 2. How currents on a surface impose a discontinuity. 
Radiating electric currents create a magnetic field which curls 
around the current thus imposing a discontinuity in the mag- 
netic field on either side of the boundary. Likewise, magnetic 
currents create an electric field which curls around the current 
imposing a discontinuity in the magnetic field 



On side one of the boundary we have some distri- 
bution of electromagnetic fields which we refer to as 
{Ei, Hi}. On side two of the boundary we have a 
different distribution of electromagnetic fields which 
are referred to as {E 2 ,H 2 }. Because the tangential 
fields are different on either side of the boundary we 
must place something at the boundary to satisfy the 
discontinuity since electromagnetic fields are naturally 
continuous. From electromagnetic theory we find that 
what is required at the interface are electric and mag- 
netic surface currents [2] equal to the discontinuity of 
the tangential fields. This is expressed as 



fix [E 2 -Ei] = -M 8 , 
fix [H 2 -Hi]=J s . 



(1) 
(2) 



The presence of these surface currents enforce the 
discontinuity in the fields as shown in Fig. 2. Without 
the currents, the fields are naturally continuous and will 
not change their propagation characteristics. However, 
by placing these currents a discontinuity can be created 

We now make a couple of observations about this 
scenario. The first observation is the active nature of 
the boundary. By this we mean that these currents are 
impressed at the boundary and radiate a field to create 
the overall discontinuity in the field. As we will show 
below for specific examples, we can envision passive 
scenarios by constraining the amplitudes of the field on 
either side of the boundary. 

Another observation is that we require both an 



electric and a magnetic current at the boundary. This is 
because we need to alter both the electric and magnetic 
field across the interface to alter the propagation of 
the wave and thus we need two degrees of freedom, J s 
and M s . This can be thought of as a Huygens source 
establishing the required discontinuity in the wave across 
the interface as it creates the desired wavefront. Also it 
is interesting to note that the discontinuity in the fields 
is enforced along a boundary only and does not require 
any bulk materials. A potentially beneficial result that 
will lead to thin conformal designs that we will further 
develop below. 

Finally, it can be noted that this idea of impress- 
ing electric and magnetic currents across an interface is 
not necessarily 'new' as it has been used in theoretical 
and numerical electromagnetics when applying the 
equivalence principle [2]. However, what is different 
here is the interpretation of these currents as physical 
quantities which can be used to engineer the propagation 
of the electromagnetic wave across the boundary. This 
then raises another questions: how can these currents 
be physically realized? To answer this we will introduce 
two ways to impress electric and magnetic currents 
along a boundary. Firstly, using active electric and 
magnetic dipoles. Secondly using passive impedance and 
admittance surfaces. 



Electric and Magnetic Dipoles 

The simplest way to implement a sheet of surface 
current, whether electric or magnetic, is to discretize 
the current using dipoles. This is a well known idea 
from antenna array theory where discrete arrays are 
approximated by a surface current [3]. Here a set of 
electric and magnetic dipole arrays is needed. Such a 
configuration has been investigated for improving the 
receive characteristics of an antenna array [4], though 
here we are looking at such combined arrays as a 
source for discontinuous wavefronts. For the electric 
currents, the current is discretized by an array of electric 
dipole antennas. For the magnetic currents, magnetic 
dipoles can be used to discretize the currents. Because 
of the assumed polarization of the fields, the electric 
currents are along the z-axis, J s = J s z and the magnetic 
currents are in the plane and tangential to the boundary 
M s = M s t, where t is the unit tangent vector of the 
boundary. Because of this, the electric dipole antennas 
lie along the z-axis and the magnetic dipole antennas 
along the t direction. 

The next question then is how many electric and 
magnetic dipoles are needed to appropriately discretize 
the continuous currents in Eqs. 1 and 2. Assuming we 
can describe our boundary as a parameterized curve 



of coordinate u , let us divide up the boundary into 
segments of length d u with a unit height d z = 1. 
Then the electric and magnetic dipole moments are 
given by p e = jujJ s d u d z , p m = jujM s d u d z . These 
dipole moments can be approximated as point sources 
spaced every d u which gives, p e , n = P&S(u — nd u ), 
Pe,n = Pe#(^ _ nd u ). The discretized electric and 
magnetic currents then are then given by a sum of these 
discrete dipole moments, 

oo 

J s ,d = JZ Pe(nd u )6(u-nd u ), (3) 

n— — co 

oo 

M sd = ^2 Pm(nd u )5(u-nd u ). (4) 

n= — oo 

The spacing of these dipole moments, d u , must sample 
the currents so that they capture the 'fastest' spatial 
variation of the currents at the boundary. Thus from 
sampling theory [5], we set d u to be small enough such 
that all the spatial variation of the electric and magnetic 
currents is captured. This can be found by taking the 
Fourier transform of J s and M s . In some cases, d u may 
be different for the electric and magnetic currents lead- 
ing to a different number of electric and magnetic dipoles. 

To physically implement these electric and mag- 
netic dipoles, the simplest way is to use small antennas. 
For the electric dipole, a traditional metallic strip of 
current can be used, while for the magnetic dipole a 
metallic loop of current would suffice. 



Impedance and Admittance Surfaces 

Another way to introduce a discontinuity in the fields 
is to use a surface which scatters the field. Such sur- 
faces can be characterized by an impedance [6]. These 
impedance surfaces form a boundary condition which 
have the following properties: 



we can envision another boundary condition which we 
will refer to as an admittance surface which is described 
by the following boundary condition: 



n x [E 2 - Ei] = 0, 
nxlHa-HxHJs, 
h x Ei = Z S J S 



(5) 
(6) 

(7) 



where in general Z s is a complex impedance that relates 
the discontinuity in the magnetic field to the continuous 
electric field. Physically these surfaces are made up of 
planar metallic scatterers on which the induced currents 
radiate a scattered field. This scattered field when added 
with the incident field yields a total field which satisfies 
Eq. 5- Eq. 7. 

Of course, given the discussion above, this bound- 
ary condition alone would not be able to introduce a 
discontinuity in the field because we need both electric 
and magnetic currents. Thus, arguing from duality [2], 



n 



Eo — Ei 



-M s 



n x [H 2 - Hi] = 0, 

n x Hi = y s M s 



(8) 

(9) 

(10) 



where Y s is a complex admittance that relates the 
discontinuity in the electric field to the continuous 
magnetic field. It is worth asking what such a surface 
would look like. In this case the scatterers would be 
metallic loop-like objects on which induced loops of 
current radiate a scattered field equivalent to what the 
magnetic current in Eq. 8 would radiate. 

With these two boundary conditions, a combined 
surface made up of both impedance and admittance 
boundary conditions can be used to introduce a discon- 
tinuity in the electric and magnetic field. To achieve 
this we can modify the scenario given in Fig. 1 by 
dividing up the fields on either side of the boundary. 
Here the electric and magnetic fields are divided into 
two groups. The first group is the continuous electric 
field and magnetic fields given by {E c i, E c2 } and {H c i, 
H C 2J. At the boundary the continuous fields are equal 
to each other with E c i = E C 2 and H c i = H C 2- The 
other group is the discontinuous fields which are given 
by {Edi, Ed2} and {Hdi, Hd2J- At the boundary these 
fields are not equal to each other. Now by grouping the 
continuous electric field and the discontinuous magnetic 
field, an impedance at the boundary can be defined 
by applying Eq. 5-Eq. 7. Likewise, by grouping the 
continuous magnetic field and the discontinuous electric 
field an admittance can be defined from Eq. 8-Eq. 10. 
This combined surface can then be used to introduce a 
discontinuity in the electromagnetic field. 

In general without placing any constraints on the 
fields across the interface, Z s and Y s will be complex 
and the real part of both quantities could be less than 
zero [7]. An impedance or admittance with a negative 
real part implies that the surface is active and supplying 
power to the system. However as we will show below 
we can envision specific scenarios where the real part is 
constrained to be passive and lossless. 



ENFORCING DISCONTINUITIES TO CLOAK 
AN OBJECT 

As we demonstrated in [1] an array of electric and 
magnetic dipoles can be used to cloak an object. Using 
this idea of discontinuous electromagnetic fields we 
introduce two other cloaking schemes which can hide an 
object by introducing a discontinuity in the field at the 



boundary. 



and magnetic currents are found which are given by, 



In the first example, a discontinuity is imposed on 
the boundary of a scatterer illuminated by a plane- wave. 
By carefully choosing this discontinuity we can create a 
cloak which cancels the scattered field outside the cloak 
while leaving the fields undisturbed inside the cloak. 
This creates a cloak which can interact with the field 
inside the volume without disturbing the field outside, 
similar to the anti-cloak proposed in [8, 9] but without 
using bulk materials. 

Our object to be cloaked is a dielectric cylindri- 



E s =0,H s =0 




-J s ,-M s ^ 



FIG. 3. A schematic of the active cloak which is constructed 
by enforcing electric and magnetic dipoles on the boundary of 
the dielectric cylinder. In this section two cloaking schemes 
are demonstrated by altering the weights of the electric and 
magnetic dipoles on the boundary. In the first example the 
fields are canceled outside of the scatterer without disturbing 
the interior fields {Ei n t,Hi n t}. In the second example the 
electric and magnetic dipoles create a scattered field exterior 
to the cloak that looks like a metallic cylinder. 



cal scatterer situated at the origin and of infinite extent 
along the z-axis. The radius of the scatterer is a = 0.7A 
at an operating frequency of 1.5 GHz while the dielectric 
constant of the material is e r = 10. The incident field 
is a plane-wave described by £^e _j/cx z. This is depicted 
in Fig. 3. The total fields outside the scatterer can be 
decomposed into an incident field and a scattered field, 
while the fields inside the scatterer are left in terms of 
the total field. At the boundary of the scatterer these 
fields are continuous. However we want to impose a 
discontinuity in the field at the boundary of the scatterer 
by adding the negative of the scattered field outside the 
object without disturbing the fields inside the scatterer. 
This gives a discontinuity that is described in terms 
of the negative of the scattered field, — E s ,— H s . By 
plugging the scattered field into Eq. 1-2, a set of electric 



M s = -n x (-E s ) = 

n—oo 

-k 2 J2 A sn H^(kp)e-^ + ^4>, (11) 

n= — oo 

J s = n x (-H s ) = 

n=oo 

jcoek Y, A s „ffi 2) '(Me- j " W+f) 2, (12) 



where A sn are the coefficients of the scattered field 
expansion in terms of cylindrical wave functions and are 
given in [2]. By imposing these electric and magnetic 
currents at the boundary of the object, the negative of 
the scattered field is radiated creating a discontinuous 
field at the boundary which leaves no scattered field out- 
side while the fields inside the scatterer are undisturbed. 
Following Section. , the currents are implemented using 
electric and magnetic dipoles. To sufficiently sample the 
currents the Fourier transform of M s and J s is taken 
and it is found that 24 electric and magnetic dipoles 
placed along the boundary of the of the scatterer are 
sufficient. Generally, the larger the scatterer, the greater 
the number of dipoles that are required to sufficiently 
sample the electric and magnetic currents. 

We can simulate this array of electric and mag- 
netic dipoles surrounding a dielectric cylinder using 
COMSOL multiphysics. The results are shown in Fig. 
where the total field of the bare dielectric cylinder is 
shown along with the total field for the cylinder sur- 
rounded by electric and magnetic dipoles. Also plotted 
is the 2D bistatic radar cross section of the cylinder 

IF I 2 

defined in two-dimensions as cfrcs — ^ i7rr \E?\ [10]- We 
can see that the scattering off of the cylinder is reduced 
while the fields inside the dielectric are undisturbed as 
shown in the inset. Such a scenario would allow for 
interesting applications in sensing where any disturbance 
in the field could be minimized while still being detected. 

The second example of a discontinuity imposed along 
the boundary of a scatterer further extends this idea 
of canceling the scattered fields without disturbing the 
fields inside. Here, along with canceling the scattered 
fields of the object without disturbing the fields inside, 
we further superimpose a set of fields outside the scat- 
terer which are the scattered fields of a different object. 
This makes the dielectric scatterer look like a different 
object. In this specific example we make the dielectric 
cylinder look like a perfectly conducting cylinder. To 
achieve this we must impose a discontinuity in the field 
through a set of electric and magnetic currents which 






FIG. 4. An active cloak which cancels the scattered field without disturbing the fields inside the dielectric. On the left is a 
plot of the total out-of-plane electric field, E z for the bare dielectric cylinder without a cloak. The inset shows the fields inside 
the cylinder. The dielectric boundary is marked with a black circle. The middle plot shows the cylinder with a cloak made up 
of electric and magnetic dipoles. We can see that the field pattern now resembles the incident plane- wave only. Note also that 
the fields in the dielectric are relatively undisturbed as shown in the inset. Finally the 2D radar cross section (RCS) is shown 
on the right indicating a decrease in the scattering off of the cylinder. 



are given by , 



M s 



-nx (-E s + E sm ) 



k 2 J2 A m H^(kp)e 



jnO+f ) 



+ 



k 2 J2 A smn H^(kp)e-^^^]4>, (13) 

n— — oo 

J s = nx (-H s + H sm ) = 

n=co 

jcvek[ J2 A m HW '(fcp)e- J ' n( * +5) - 

n= — oo 
n=co 

J2 A smn H^' (kp)e-^ + ^]z, (14) 

n= — co 

where E sm and H sm are the scattered fields of a 
perfectly conducting cylinder of the same radius as 
the dielectric scatterer and A srnn are the scattering 
coefficients of the perfectly conducting cylinder, given 
by A smn = -J n (ka)/(k 2 H 2 (ka)) [2]. These electric 
and magnetic currents create a field outside of the 
cylinder which looks like the field of a perfectly con- 
ducting cylinder while canceling the scattered field 
of the dielectric object along without disturbing the 
fields inside. To implement the currents, electric and 
magnetic dipoles are used again and like the pre- 
vious example, 24 dipoles of either kind are required 
to sufficiently sample the currents for the same geometry. 

This configuration is also simulated using COM- 
SOL multiphysics and the results are shown in Fig. 5. 
Here the total field of a bare perfectly conducting 
cylinder is shown along with the total field of the 
active cloak. Here we can see from the electric field 
pattern, we have made a dielectric cylinder look like a 
conducting cylinder. Again the field inside the dielectric 
cylinder is undisturbed. Finally the 2D bistatic radar 
cross sections is shown and we can demonstrate a good 
agreement between the scattering of a metal cylinder 
and our cloak and the noticeable difference between a 



bare dielectric cylinder. This demonstrates the ability of 
this array of electric and magnetic dipoles to disguise a 
dielectric cylinder as a metallic cylinder by imposing a 
discontinuity in the field at the boundary. 



Comments 

With these basic examples it becomes clear that 
being able to impose a discontinuity in the electro- 
magnetic field is a powerful concept which allows for 
the realization of some novel applications. We can see 
that by imposing electric and magnetic currents at a 
boundary through discrete dipole arrays, the wavefronts 
can be manipulated in many possible ways. As we 
have demonstrated so far by having both electric and 
magnetic dipoles we are able to radiate a wavefront 
which creates our desired field as per Huygens principle. 

A distinguishing feature of this proposed approach 
stems from the need to know the scattered field and 
thus the incident field a priori. This is inherent in this 
approach of creating discontinuous fields, as the discon- 
tinuity must be inserted into a known field distribution 
[11]. However for applications where the incident field is 
known or can be determined using other methods, this 
approach shows significant promise. 



ENFORCING DISCONTINUITIES TO REFRACT 
AND REFLECT A PLANE WAVE 

Another interesting example of wavefront manipula- 
tion is through altering the direction of a plane wave 
either through refraction or reflection. As stated earlier, 
this is typically done with bulk materials which alter the 
fields by forcing continuity across the boundaries. Here 
we now apply the idea of creating a discontinuity in the 
electromagnetic field to refract a plane wave. 

The basic scenario is described in Fig. 6 where a 






FIG. 5. An active cloak which cancels the scattered field of the dielectric cylinder while making the dielectric cylinder look like 
a metallic cylinder. On the left is a plot of the total out-of-plane electric field, E z for a bare metallic cylinder. The boundary is 
marked with a black circle and is the same size as our dielectric cylinder. The middle plot shows the dielectric cylinder with a 
cloak made up of electric and magnetic dipoles. We can see that the field pattern now resembles the total field of our metallic 
cylinder, which disguises the dielectric cylinder. Note also that the fields in the dielectric are undisturbed as well. Finally 
the 2D radar cross section is shown on the right demonstrating quantitatively how the dielectric cylinder resembles a metallic 
cylinder when cloaked with the active dipole array. 




FIG. 6. Refraction and reflection at an interface. If the in- 
cident plane-wave is refracted or reflected at the boundary 
in free-space, a discontinuity exists at the interface which re- 
quires an electric and magnetic current at the interface. We 
can impose those currents using a discrete active electric and 
magnetic dipole array or we can use a passive impedance and 
admittance surface, provided we constrain the amplitudes of 
the fields. 



plane wave in free-space is incident on a boundary 
located at x — 0. At this boundary, a discontinuity in 
the field is imposed and the wave suddenly bends as 
it traverses the boundary. To design the boundary we 
will look at two ways to impose the discontinuity using 
active electric and magnetic dipoles as well as a passive 
impedance/admittance surface. 



Refraction Using Active Arrays of Electric And 
Magnetic Dipoles 

To use an array of electric and magnetic dipoles to 
refract a plane-wave wave we can simply subtract the 
differences between the incident and refracted plane- wave 
at the surface. Using Eqns. 1 and 2 we can find the 
electric and magnetic currents which are given by 

M s = -y [E t e- jkxcos9t - E ie - jkxcos ^} , (15) 

-cosO t E t e- jkxcos6t --coseiEie- jkxcos h^) 
i) V J 



where E{ and 0i are the amplitude and direction of the 
incident wave and E t and t are the amplitude and 
direction of the refracted wave and ri = x — is the 

impedance of free-space. Here we have assumed that 
there is no reflected wave generated by the interface. 
We also note a similar solution could be found for a 
reflected wave only. Following again the procedure given 
in Section , an array of electric and magnetic dipoles 
can be constructed at the x = plane. Here we need 
to place an electric and magnetic dipole every 2A/3 to 
sufficiently sample the continuous electric and magnetic 
currents. 

Again we can simulate this scenario using COM- 




FIG. 7. An array of electric and magnetic dipoles at a surface 
which interfere to generate a negatively refracted plane- wave. 
The total electric field along the vertical axis is plotted. The 
incident field is along the 0i = 20° direction and the refracted 
wave along the 0t — —20°. 

SOL multiphysics. In our specific example we generate 
negative refraction where a plane-wave incident at an 
angle of 6i = 20° refracts to an angle of 6 t = —20°. 
This is depicted in Fig. 7. This demonstrates the basic 
concept of a discontinuous electromagnetic field at a 
planar interface. Again we can see how both electric and 
magnetic currents are required to radiate the refracted 
wavefront. 

For refraction, it is worth asking if this is possible 
using a passive method like in 'normal' refraction. Thus 



we will investigate the use of impedance and admittance 
surfaces to refract an incident plane- wave. To ensure the 
passivity of the impedances we will look at constraining 
the amplitudes of the refracted fields. 



Other Ways Of Refracting A Plane Wave At A 
Surface 



It is worth pausing here to discuss other methods 
of refracting a plane wave using a surface as this is a 
very active area in the literature. In the microwave 
regime, transmitarrays and reflectarrays have been used 
as scattering surfaces to bend or reflect incident waves 
into arbitrary directions as they pass through a surface 
[12, 13]. These surfaces can often be characterized as 
impedance surfaces [13]. For these microwave devices 
multiple surfaces are often used to optimize bandwidth 
and/or to minimize reflections at the surface, however 
the main concept can be reduced to a single surface. 
At the infrared and optical wavelengths similar designs 
have been constructed using either passive arrays of 
nano-antennas [14, 15] or blazed grating-like structures 
[16]- 

Despite the different language used between the 
microwave and optical communities, these ideas all 
follow a similar logic that is easily described using 
a refraction law derived from a phase discontinuity 
between the refracted and incident fields [14, 17]. Here 
a linear phase shift along the elements that make up the 
surface are used to alter the direction of the incident 
plane- wave. This is expressed as 



smt/ t 



sin Uj = 



k dy' 



(17) 



where $ = ky {sin Ot — sin Oi) is the phase difference 
along the surface between the incident and refracted 
wave. Thus if we have a linear phase variation along 
a surface, we expect an incident plane- wave to refract 
along the t direction after passing through the surface. 
This idea follows somewhat analogously from linear 
antenna array theory [18]. However, as demonstrated 
in [13], the problem with structures that provide a 
linear phase-shift along a surface is that they always 
excite other modes (other plane-waves or Floquet modes 
for periodic structures) which propagate in different 
directions other then the desired angle of refraction 
t . And if one examines the simulated or measured 
results in the papers listed above, this behaviour can be 
observed. 

There are two changes which can be made to these 
designs to get around this problem. The first and most 
fundamental change is to impose both electric and 
magnetic currents to create a proper discontinuity in 



the electromagnetic field as we have proposed. The 
second change is that the elements which make up 
our scattering surfaces do not vary linearly to impose 
a linear phase shift in the electromagnetic field. We 
will now demonstrate this using both impedance and 
admittance surfaces to refract a plane- wave. 



Refraction Using Passive Impedance And 
Admittance Surfaces 



As stated, impedance and admittance surfaces can be 
used to implement a discontinuity in the electromagnetic 
field. In line with the previous work described in the 
last section we would like to make our impedance 
and admittance surfaces passive and lossless, implying 
5R(Z S ) = and Sft(Ys) = 0. This then requires us to 
carefully construct the amplitudes of the waves which 
scatter off of the surface. To find the required Z s and 
Y s to refract a plane- wave, we will break up the fields 
into a set of continuous electric fields and discontin- 
uous magnetic fields and vice versa as stated in Section. . 

For both the impedance and admittance surfaces 
illuminated by a plane-wave we can describe a set of 
fields which scatter off of the surface and satisfy the 
boundary conditions given in Eq. 5-Eq. 10. Because of 
the symmetry across the y-axis, both the impedance 
and admittance surfaces scatter the incident wave into 
plane-waves which propagate from either side of the 
boundary in the same direction. This forces us to 
analyze the problem by constructing a symmetric set of 
scattered fields on either side of the surface. 
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FIG. 8. A plane- wave incident on an impedance surface. 
Shown are the incident and the scattered fields. On side 2 of 
the boundary the total fields are the scattered fields summed 
with the incident plane wave (not-shown). 



Recalling that our ultimate goal is to arrive at a 
total field which refracts the incident plane-wave across 
the boundary, our combined impedance and admittance 
surface has to do two things. First it must cancel the 
incident field that exists beyond the interface on side 
two. Second it must scatter the field into the refracted 
beam on side two as well. Thus an incident plane-wave 



k H t| 



Ei Ej sm r^ k "^ :e k 



H-r^K \ H ism 

1 :«\ * 



z tm 



Z y © 

I I I I I I I I I I I 1:1 + II H V I I I ,1 II II 

Hi«m ! R.__ | Y s © 



'ism 



-ism 



m ^ 

■ i K v 



c tm ▼ 



:e t ^k 



""tm ' 



FIG. 9. A plane- wave incident on an admittance surface. 
Shown are the incident and the scattered fields. On side 2 of 
the boundary the total fields are the scattered fields summed 
with the incident plane wave (not-shown) . To create a surface 
which primarily reflects the incident plane-wave as opposed 
to refracting it, the sign of E tm should be flipped. 
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scattering off of either the impedance or admittance 
surface must generate four-scattered plane- waves, two 
on either side of the interface. This is pictured in 
Fig. 8 for the impedance surface and in Fig. 9 for the 
admittance surface. These scattered fields in Fig. 8 
form a set of continuous electric fields and discontinuous 
magnetic fields which will define the impedance surface. 
The scattered fields in Fig. 9 form a set of continuous 
magnetic fields and discontinuous electric fields which 
correspondingly define the admittance surface. 

We can now use these fields to analyze the impedance 
and admittance boundary condition respectively. Start- 
ing with the impedance boundary condition we plug 
these fields shown in Fig. 8 into Eq. 7 to give us 
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We can now solve for E ise and E te , which are the ampli- 
tudes of the scattered fields, by forcing $l(Z s ) = 0. This 
will also give us an expression for $s(Z s ) = X s . Doing 
this we find that 
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Ej cos 0j 



cos 6i + cos 6 t 
And our expression for the surface reactance is, 
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where $ = ky (sin 8 1 — sin^), the phase shift between 
the refracted and incident plane waves. Note that our 
expression for X s is not linear with respect to y (or 
<l>), implying that we need a non-linear gradient in 
our impedance surface to generate a linear phase-shift 
for our refracted wave. Also note that this is the best 
possible case for a single impedance surface as we can 
channel the incident plane- wave into a minimum of four 
other plane- waves, and thus only succeed in partially 
refracting the beam. 

We now do the same for the admittance surface at 
x = 0. Taking the fields in Fig. 9 and inserting them 
into Eq. 10 and forcing Sft(Ys) = we get 
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Again we also get an expression for the surface suscep- 
tance, 
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which is a non-linear function of y, the spatial coordinate 
along the boundary. Note also that both Eq. 21 and 
Eq. 23 depend on <£, the phase shift between the incident 
and refracted wave. 

As we have stated before, the complete effect 
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FIG. 10. A plane- wave incident on a combined impedance 
and admittance surface. Shown are the incident and the total 
fields on either side of the boundary. Because the impedance 
and admittance surface create induced electric and magnetic 
currents we can successfully refract the incident plane-wave 
without any other fields on side two of the boundary. Also 
note that the majority of the power is scattered into Et as 
the amplitude of Ei S is much smaller as shown in Fig. 11 

comes when our impedance and admittance surfaces 



are superimposed creating a discontinuity in both the 
electric and magnetic field. Superimposing the scattered 
fields off of both the impedance and admittance surface 
we get the following set of fields shown in Fig. 10. We 
now have an incident plane-wave which when scattering 
off of both admittance and impedance surfaces gives us 
a refracted field with an amplitude given by, 

^ ,- ^ 2Ei cos 0i . A N 

E t = E te + E tm = * \ . 24 

cos Oi + cos t 

We also have a reflected field which reflects off the surface 
in a specular fashion with an amplitude of, 

Ei(cos9i — cos t ) 



cos Oi + cos t 



(25) 



Note however, that this reflected field is much smaller 
than the refracted field. 

To further show that these are the only fields gen- 
erated by combined impedance/admittance surface, we 
can look at the power entering and leaving the scattered 
surfaces by integrating the fields over an area enclosed 
by the surface as shown in Fig. 10. The power entering 
and leaving the surface is given by, 



P=^ff»[ExH*]- MS. 



(26) 



For plane-waves this reduces to the power in the normal 
component of the field that enters or leaves a surface with 
area, A, which gives, 
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\Ei\ 2 cos Oi - E 2 t cos t - Ef s cos 0*) = (27) 



which demonstrates that all the power from the incident 
plane-wave is scattered into the refracted and reflected 
beam. 

An equivalent setup for an impedance and admit- 
tance surface that reflects an incident plane- wave along 
an arbitrary angle can easily be extended from the 
derivation given here by simply negating the direction of 
E tm in Fig. 9 and re-deriving the subsequent equations. 

We can summarize these results in a set of design 
curves which describe how a plane- wave refracts through 
a combined impedance and admittance surface. This 
is plotted in Fig. 11 and describes the solutions of 
E t , Ei S , X s and B s for all possible combinations of 
Oi and t . The plots on the top- left and top-right of 
Fig. 11 give us the amplitude of E t and Ei S respectively 
for different values of t and Oi. Note that for most 
reasonable angles of incidence the amplitude of the 
refracted beam, E t is close to 1 while the reflected 
beam is close to showing that our combined surface 
is ideal for refraction with minimal reflection. In the 
bottom-left and bottom-right plots we see contours of 





FIG. 11. Design curves summarizing refraction through an 
impedance and admittance surface for an incident plane wave 
of amplitude E% = 1. All plots are drawn with incident angle, 
Oi, on the horizontal axis and the refracted angle, t , on the 
vertical axis. The top-left plot gives us the amplitude of the 
refracted beam, E t while the top-right plot gives us the am- 
plitude of the reflected beam Ei S . The bottom two curves are 
contour plots of the impedance and admittance of the screen 
at a given point, (y = A) for different incident and reflected 
angles. 



constant impedance and susceptance respectively for 
different combinations of t and Oi. Note here, that for 
small values of both t and Oi the curves overlap. This 
can be seen from inserting small-angle expressions for 
the angles into Eq. 21 and Eq. 23. This implies that for 
small angles, a fixed combined impedance/admittance 
surface can predictably refract an incident plane wave of 
an arbitrary angle of incidence. However, as both these 
values diverge the impedance and admittance curves 
diverge from each other as well. In general then, like the 
active cloaks shown here and in [1], we are constrained 
with knowing the incident field to design our surface for 
a desired angle of refraction, Of. 

Again, these admittance and impedance surfaces 
act like a passive Huygens source which generates the 
required wavefronts to construct the refracted field. 
The incident plane-wave along a direction Oi induces 
electric and magnetic currents on the impedance and 
admittance surface respectively which independently 
radiate a scattered field which interferes to form a 
refracted plane- wave in the direction of t 

If we had not been so careful as to constrain the 
amplitudes of the refracted fields we would have arrived 
at impedance and admittance surfaces which had a 
non-zero and possible negative real part. This would 
have implied a surface which scattered more power 
then was incident on it and would be analogous to the 
example in Section where we used an active dipole 
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array to refract a plane- wave. As we have shown for 
refraction, the power in the incident field can simply be 
redirected into the refracted field by a passive surface as 
long as the amplitudes are constrained. Thus an active 
surface is not truly needed (though it does demonstrate 
the concept of a field discontinuity). To contrast, for 
cloaking the power in the incident and scattered field is 
already conserved [19]. Thus to enforce a discontinuity 
in the field which cancels the scattered field the power 
must come from elsewhere, in this case an array of 
electric and magnetic dipoles (Note that this is different 
than bending the light around the object as in [20] or 
resonating out a specific multipole order of the scattered 
field [21] which are passive cloaking schemes but which 
do not superimpose a field into the problem). 

Finally, we note that this construction of impedance 
and admittance surfaces gives the full description of 
how to refract a plane- wave across a surface. As a 
point of comparison, we can look at how refraction and 
reflection of a plane-wave across a material interface 
is derived. Here we notice three attributes about this 
process. First it is derived from the boundary conditions 
at a material interface. Secondly, it tells us how the 
field will refract and reflect (Snell's Law). Third, it tells 
us the amplitudes of the field (the Fresnel reflection 
coefficients). Here we have the same attributes but re- 
alized using a combined impedance/admittance surface, 
for a surface. Starting from the boundary conditions 
given by impedance and admittance surfaces, we have 
derived how the field will refract across a surface for 
different impedances (Eq. 21 and Eq. 23) and we have 
found the amplitudes of those fields in Eq. 24 and 
Eq. 25. This gives us a complete picture for how to tai- 
lor a surface to refract (or reflect) an incident plane- wave. 

In comparison, the generalized refraction law in 
Eq. 17 simply tells us the phase gradient between the 
two refracted plane- waves. It does not tell us how to 
synthesize a surface. (And as shown, trying to synthesize 
a surface which simply mimics the phase gradient is 
incomplete). Thus this equation is descriptive, as it 
describes how a plane-wave refracts across a surface, but 
it is not prescriptive in that it does not tell us how to 
design our surface. Eq. 17 is useful as a preliminary step 
as we can envision an arbitrary gradient along a surface 
which encodes some functionality (focusing, beam 
steering, vortexes) with respect to some incident field. 
However we would then need to use the equations given 
in Eq. 21 and Eq. 23 to synthesize the required surface 
using both electric and magnetic currents generated by 
the impedance and admittance surfaces. 



Examples 

We can verify these results by using two simple full- 
wave simulations. The first example is to use the two- 
dimensional method of moments to solve for the required 
impedance and admittance as well as the fields radiated 
by the surface. The second example is using a commer- 
cial full- wave solver (HFSS) to demonstrate a very basic 
physical implementation of this concept above. 



Method of Moments Verification 

To analyze the refraction through a surface we can 
use a two-dimensional method of moments procedure 
to numerically synthesize the required admittance and 
impedance as well as to find the radiation from the 
induced currents on the surface, confirming the results 
above [7, 10]. 

For the impedance surface the induced electric cur- 
rents on the surface can be found from examining the 
relationship between the incident, scattered and total 
fields at the boundary 



EtotU=o — Ei| x =o + E s | x =o- 



(28) 



The scattered field is created by the radiation of the elec- 
tric currents induced on the surface which is given by [10], 



UJ\1 



r Z J s (y')H 2 (k\y-y'\)dy\ (29) 

J-L/2 



r L/2 



where L is the length of the impedance surface. 

Likewise, for the admittance surface, we can solve 
for the magnetic currents by examining the incident, 
total and scattered field at the boundary 



nx Htotl 



nxHil 



h x H s 



(30) 



Here we must find the tangential magnetic field created 
by the radiation of the magnetic currents on the surface 
which is given to be, 



Ei, 

V 



n x H to tU=o = — cos^e-^ Bin ^x 



\uj\± \ dy 2 



f M s (y')H 2 (k\y-y'\)dy%31) 

J-L/2 



Both Eq. 29 and Eq. 31 can be solved for J s and M s 
respectively using the method of moments since we know 
both the desired total fields and incident field at the 
boundary. To solve, the impedance surface is discretized 
into N segments and is solved for using pulse-basis 
functions and point matching [10]. Once J s and M s are 
found, these can be respectively inserted into Eq. 7 and 
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screen we have a field which resembles a plane- wave prop- 
agating along 6 t = 30° . On side one of the screen there 
is a small reflected field as predicted by Eq. 25 which is 
plane-wave along the 8{ = 0° direction. Note however 
that the amplitude of these fields is much smaller than 
the refracted fields on the other side of the screen. This 
indicates that a surface of induced magnetic and electric 
currents through impedances and admittances allows for 
a discontinuity in the electric and magnetic field to be 
realized. 



Physical Implementation 



FIG. 12. Method of moments verification of an impedance 
and admittance surface to refract an incident plane- wave. 
Solid, coloured curves are calculated method of moments re- 
sults. Dashed black curves are theoretical results from the 
previous section. The top-left and top-right figures are the 
calculated reactance and susceptance of the two surfaces. We 
can see good agreement with the theoretical results as well as 
their non-linear dependence on the spatial coordinate of the 
surface. The bottom- left and bottom-right plots show the cal- 
culated and theoretical electric field on side-two and side-one 
of the surface respectively (real and imaginary parts in red 
and blue respectively ). Note the much larger amplitude of 
the the refracted beam compared to the reflected beam. 



As stated before, the impedance and admittance 
surfaces set up induced electric and magnetic currents 
to impose a discontinuity in the field. To physically 
implement such a surface we can take a cue from the 
active version of this idea and use passive dipoles and 
loops to create our induced electric and magnetic cur- 
rents. We will demonstrate this idea by using an array 
of loops loaded with a reactive impedance to create our 
induced magnetic current. We will use dipoles loaded 
reactively to create induced electric currents. By varying 
the reactive loading on each loop/dipole the impedance 
and admittance of the surface can be tuned. 



Eq.10. 



To demonstrate this we take an example where our 
incident field is in the normal direction, 6{ = 0° and 
our desired refracted beam is along t = 30°. Our 
impedance and admittance surfaces are L = 10A long 
at / = 1.5 GHz. We divide our surface up into N = 1000 
segments that are A/ 100 long. The calculated reactance 
and susceptance are plotted in Fig. 12 along with the the- 
oretical values given by Eq. 21 and Eq. 23 where good 
agreement can be seen between the two sets of curves. We 
can also find the total field on either side of the screen by 
finding the fields radiated from both sets of currents. On 
side two of the boundary these fields must be summed 
with the incident field. These radiated fields are given 

by, 
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Plotting these radiated fields in Fig. 12, we find a field 
that resembles the theoretical prediction fairly well on 
either side of the screen with the discrepancies coming 
from the finite nature of the screen. On side two of the 



To model this we will use Ansys's HFSS. Again we 
will use the same example as before at 1.5 GHz with 
an incident field at 6i = 0° and a refracted field at 
t = 30° . We will place our simulation in a parallel-plate 
waveguide to simplify the computation and emulate a 
2-D environment. Taking the impedances and admit- 
tances given in Fig. 12, each loop/dipole is designed to 
implement the desired impedance at its location along 
the surface, with the loops/dipoles spaced every A/10. 
This is shown in Fig. 13. To illuminate the surface a 
Gaussian beam is used with a 3A focal spot designed 
to occur at the surface. With this design we offer a 
couple of caveats. First the unit cells are not optimal 
and better designs are possible. Secondly the impedance 
variation along the surface is not optimized for a finite 
beam such as a Gaussian beam. Nonetheless the main 
point can be demonstrated here with this simple array 
of passive loops and dipoles. 

[6Z) rj^ e s i mu i a ted results are shown in Fig. 13 where 
the field-plot at 1.5GHz is shown. Here we can see that 

gagthe Gaussian beam is refracted by the surface from an 
incident direction of 8i = 0° to a refracted direction 

(34<^f 6 t = 30°. Reflections are also minimal though any 
discrepancies can be attributed to the reasons given 
above. The far-field pattern of the scattered field is also 
plotted in Fig. 13 where we can see the main refracted 
beam at = 30°. The peak at = 0° is the field radiated 
by the surface to cancel out the incident field on side 
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two of the surface (as we are plotting the scattered 
field). There is a reflected peak at 6 = 180° due to the 
small reflected field that is generated by the surface 
and a reflected peak at = 150° that is not completely 
canceled out due to the imperfections in the surface. 
This can be improved upon further optimization of the 
screen. 



CONCLUSION 

Through the examples presented in this work, we 
have demonstrated how the idea of discontinuous fields 
can be used to alter, shape and control electromagnetic 
wavefronts. By placing electric and magnetic currents 
at a boundary, the fields can be altered in a variety of 
novel ways including the ability to cloak an object and 
the ability to refract (or reflect) an incident field. We 
have shown that we can create our electric and magnetic 
currents in one of two main ways, either by placing a 
discrete array of active electric and magnetic dipoles 
which impress a discontinuity in the field. Or by using 
a passive impedance and admittance surface, on which 
electric and magnetic currents are induced to create a 
discontinuity. 

This approach of imposing a discontinuity in the 
electromagnetic field opens up many new possibilities 
across the electromagnetic spectrum. Some interest- 
ing areas include the synthesis of active electric and 
magnetic dipole arrays as well as passive impedance 
and admittance surfaces at both microwave frequencies, 
using well known technologies in phased-arrays and 
impedance surface [3, 22], and at optical and THz 
frequencies using more recent developments such as 
nano-antennas, nanophotonic phased-ararys and nano- 
circuits [23-25] . The design of novel devices using 
discontinuous fields including lenses, cloaks, reflectors, 
can also be envisioned using this new technique and 
presents a promising avenue for future work. 
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FIG. 13. Simulation results from HFSS for a physically realized impedance and admittance screen. On the left a schematic of 
the surface. The impedance screen is made of react ively loaded dipoles with the reactive loading shown in red. The admittance 
screen is made of react ively loaded loops with the loading shown in blue. In the middle is a plot of the total vertical electric 
field for a simulation of the impedance screen for an incident Gaussian beam. The plot on the right is a far-field plot of the 
scattered electric field. 



